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ABSTRACT 

We are carrying out multi-frequency radio continuum observations, using the Aus- 
tralia Telescope Compact Array, to systematically search for collimated ionized jets 
towards high-mass young stellar objects (HMYSOs). Here we report observations at 
1.4, 2.4, 4.8 and 8.6 GHz, made with angular resolutions of about 7", 4", 2", and 1", 
respectively, towards six objects of a sample of 33 southern HMYSOs thought to be 
in very early stages of evolution. The objects in the sample were selected from ra- 
dio and infrared catalogs by having positive radio spectral indices and being luminous 
(Lboi > 2 x 10 4 Lq), but underluminous in radio emission compared to that expected 
from its bolometric luminosity. This criteria makes the radio sources good candidates 
for being ionized jets. As part of this systematic search, two ionized jets have been dis- 
covered: one previously published and the other reported here. The rest of the observed 
candidates correspond to three hypercompact HII regions and two ultracompact HII 
regions. The two jets discovered are associated with two of the most luminous (7 x 10 4 
and 1.0 x 10 5 L@) HMYSOs known to harbor this type of objects, showing that the 
phenomena of collimated ionized winds appears in the formation process of stars at least 
up to masses of ~ 20 M & and provides strong evidence for a disk-mediated accretion 
scenario for the formation of high-mass stars. From the incidence of jets in our sample, 
we estimate that the jet phase in high-mass protostars lasts for ~ 4 x 10 4 yr. 

Subject headings: ISM: jets and outflows — radio continuum: stars — stars: formation 
- stars: individual (G337.4032-00.4037) 



1. INTRODUCTION 



The determination of whether massive stars (M > 8M & ) are formed by an a ccretion process 
similar to that inferred for low- mass stars, with an accreting disk (jShu et al.lll987l ). or via merging 
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processes (e.g., iBonnell et al.lll998l ) is one of the main observational challenges in the field of star 
formation. Theoretically, an accretion disk phase could even be more important for high-mass 
stars, compared to the low -mass case, in term s of the fraction of the mass of the central protostar 
accreted from such a disk (jKuiper et al.ll2010l ). 



One of the striking features of the disk accretion theory is its connection with the ejection of 
a collimated jet. The consensus is that detecting the presence of such jet is a sufficient condition 
to ensure that disk accretion is taking place. Whether jets are key in al lowing disk material to lose 
angular momentum, and fall onto the protostar, is still an open issue ([Livid |2009| ) . However, the 
latter idea is suggested by theoretical work in magneto-hydrodynamic modeling of the disk accretion 
and jet ejection process. In addition, jets form a natural candidate to explain the momentum and 



energy source of the com monly observed massive bipolar molecular outflows (e.g. iBeuther et al 



2002 



Kim Kurtal2006l ). Collimated jets detected toward O-type forming stars would be strong 
evidence favoring a disk-mediated accretion process. On the other hand, if t hey are formed via 
merg ing of lower-mass stars then neither accretion disks nor jets are expected (jBally &: Zinnecker 
2005). The study of the incidence of jets in high-mass young stellar objects (HMYSOs) will then 
aid in identifying the dominating mechanism of massive star formation. 

There exists evidence that stars wit h masses up to ~ 10 (or early B-type star) are formed i n 



a disk mediated accretion scenario (e.g.. lGaray fc Lizano 



1999 



Patel et al 



2005 



Chini et al.ll200fih 



and there is gro wing observational evidence of a disk - mediated accretion c ollapse for even more 
massive objects ( Cesaroni et al. 2007 : Kraus et alJboiO ; Preibisch et al The evidence for jets 

from massive YSOs is rather scarce, however. To date there are only a handful of H MYSOs known 



to be associated with highly collimated jets and/or Herbig-Haro (HH) objects (see iGuzman et al 
2010) . Most of them have bolometric luminosities smaller than 2 x 10 4 Lq corresponding to that of a 



B0 zero-age main sequence (ZAMS) star. There is only one O-type YSO with bolometric luminosity 



Garav et al. 


2003; 


Brooks et al. 


2007) 



16547-4247, L ~ 6.2 x 10 4 L Q ; 
lack of young massive stars with spectral types earlier than B0 associated with jets and/or disks is 
an intrinsic property of the most massive stars or due to observational disadvantages. Massive stars 
are rarer and their evolutionary time scales are much shorter than those of low mass stars. With 
the currently available small sample, made from a collection of individual serendipitous studies, it 
is difficult to characterize the jet phenomena in the process of high-mass star formation. 

The presence of jets can be observationally established from high angular resolution radio 
continu um observations either directly through their elongated structures and distinctive spectral 
indices ([Reynolds! Il986l ) or indirectly by the detection of phenomena intimately associated to the 
jets, such as Herbig-Haro objects or aligned radio lobes. In this paper we report the first results of 
a systematic search for ionized jets, made using the Australia Telescope Compact Array (ATCA), 
towards a sample of HMYSOs thought to harbor this type of objects. In <J2] we describe the 
selection criteria adopted to build a sample of jet candidates selected from two catalogs of radio 
sources associated with either HMYSOs or ultracompact HII regions (UCHII). In $3] we describe 
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the ATCA observations toward six of the jet candidates, and in $3] we present the results and 
physical interpretation of the data obtained toward these six sources. 

We discuss in £|5]the search outcome and analyze the consequences of the observed jet incidence. 
In the analysis we include the results obtained toward one of the jet candidates, G345. 4938+01. 4677 
(also IRAS 1 6532—3959). A detaile d account of its characteristics and physical nature was already 



presented by iGuzman et al.l (|2010i ). Summarizing, the radio emission corresponds to a string of 
radio sources consisting of a compact, bright central component, and four aligned and symmetrically 
located lobes. The central object corresponds to a ionized jet and the emission from the lobes arises 
in shocks resulting from the interaction of the collimated wind with the surrounding medium. We 
stress that this source, which is currently one of the most massive and luminous (7.0 x 10 4 L@) 
HMYSO known to be associated with an ionized jet, was discovered as part of the systematic search 
reported in this paper. 

Finally, we summarize our work and conclusions in £j6) In the Appendix we report the results 
of observations obtained toward a sample of five radio sources which were selected in an early stage 
of this study. Three of them are radio sources associated to young high-mass stars, but they do 
not fulfill our selection criteria to be a jet candidate. 



2. SOURCE SELECTION 

In this section we describe the steps followed to build our sample of jet candidates towards 
high-mass young stellar objects. Three requirements were imposed in selecting the targets: 



Posit ive radio continuum spectral indices. First, and as expected for thermal jets ([Reynolds 



19861 ). we considered radio sources in the literature with positive spectral indices at radio 



wavelengths. 

• Association with luminous infrared sources. Since we are searching for jets associated with 
high-mass YSOs, we expect the driving sources to be luminous and still be enshrouded in 
a dense cloud of gas and dust. Thus, most of their luminosity should be re-emitted in the 
infrared and far-infrared (FIR). A second requisite to be fulfilled by our targets is that they 
be associated with luminous IRAS sources. 

• Under luminous radio objects. We expect radio jets to be only visible in a very early stage of 
evolution prior to the ultracompact H II region phase. In such an early stage the UV photons 
from the central protostar are likely t o be trapped by the infalling gas and therefore quench 
the development of an UCHII region (|Yorkelll986l : IWolfire Cassinellilll987l ). Thus the third 
requisite is that the radio luminosity of the targets be considerable lower than that predicted 
from the bolometric luminosity. 
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(2007a 


) and 


Walsh et al. 


(1998). 



rations, at 4.8 and 8.6 GHz using ATCA, toward 826 Red MSX So urces CRMS. lUrauhart et al 



2008a|). RMS sources are HMYSO candidates selected by t heir colors (ILumsden et al.ll2002f ) in the 
Midcourse Space Experiment (MSX) and 2-MASS bands (jSkrutskie et al.l 120061 ) . We considered 
all the radio sources that are within 18" (MS X beam) from the peak position of the associated 
MSX source, which amounted to 239 objects. IWalsh et al.1 ()1998l ) reported radio continuum ob- 
servations at 6.7 and 8.6 G Hz, using ATCA, towards IRAS sources with colors of UCHII regions 
(jWood fe: Churchwelllll989i ) and associated with methanol masers. They reported emission at 8.6 
GHz from 177 ultracompact HII regions. 

To be in accordance with the first requirement, we selected from the above radio sources those 
which exhibit positive radio continuum spectral index between the two observed frequencies. We 
find that approximately a 40% of the radio sources from both catalogs fulfill this first requirement. 
We considered sources with upper limits (or non-detections) at the lower frequency. In the case of 
Urquhart et al .1 ([2007a]) the sensitivity in the two observed bands is similar, hence a n on-detection 



i n the lower frequency band implies a true positive spectral index. In the case of IWalsh et al 



(j 19981 ). the observations at 6.7 GHz were made with a bandwidth of 4 MHz, much smaller than 
that used for at 8.6 GHz (128 MHz), and therefore are considerably less sensitive than those at the 
high frequency. In this case, a non detection does not necessarily implies a truly positive spectral 
index. 

To comply with the second requirement, we then selected from the above sub-sample of radio 
sources those that are associated — angular separation no larg er than 25" — with a n IRAS source. 
As distances we adopted the kinematical dis tances reported bylFaundez et alj (120041) (derived from 



observations of the CS(2— >1) transition in iBronfman et al.l 11996) or by lUrauhart et al.l (|2007bl . 



2008bl ) (derived from observations of the 13 CO(2— >4) line). In iFaundez et al.l and in the RMS 
catalog, the near-far distance ambiguity fo r sources located insid e the solar circle (R^ = 8.5 kpc) 
has been resolved for several sources (see iMottram et al.ll2011al : lUrquhart et al.M2012l ). For the 
sources in which the ambiguity remains, we adopted the near distance. Fro m the IRAS fluxes and 
the distance, we estimated the bolometric luminosity using the expression (jCasoli et al.lll986l ) 



LiRAi 



5.44 
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079 ~2~ ~33 ITsT 



(— 

Vkpc 
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where Ljras is the IRAS FIR luminosity, D is the distance and F12, F25 Fqo, and F100 correspond 
to the IRAS fluxes in Jy measured at the 12, 25, 60, and 100 fim bands, respectively. From the 
sources with positive radio spectral index, approximately 37% have molecular line observations and 
FIR luminosity Ljras > 2 x 10 4 L q . 

To comply with the third requirement, we computed the monochromatic radio luminosity, 
4ttD 2 F u , expected from an homogeneous, optically thin HII region that is excited by a single star 
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with a bolometric luminosity, Liras, using the expression (|Spitzerlll998l ) 
'4nD*F v \ ( iV* \ / v \-o-ii/ T 



Jy kpc 2 



0.12 



10 45 s -l 



?.6 GHz 



\ 0.38 

8000K J 



(2) 



where iV* is the rate of ionizing continuum photons e mitted by a ZA MS star with that luminosity 
taken from standard stellar atmospheres models (e.g.. |Panagialll973l ). In Eq. ([2]) v is the observing 
frequency, and T is the electron temperature of the ionized gas. Then we selected from the above 
subsample those radio sources for which the observed radio luminosity is smaller than that predicted 
from the latter equation by a factor of at least 10, assuming T = 8000 K. Note that Eq. ([2]) implicitly 
assumes that all the ionizing photons are absorbed by the plasma. This third and last requisite 
reduces the number of jet candidate HMYSOs approximately in 50%. 

While we argue that the third requirement is needed to select ionized jets, there are other 
explanations for the radio under-luminosity: (i) The radio emission at the observed frequency may 
arise from an optically thick H II region, in which case the radio luminosity is not proportional to 
the number of ionizing photons; (ii) Dust within the H II region absorbs an important fraction of the 
ionizing photons, (iii) The distance to the source is overestimated, (iv) The FIR luminosity is due to 
multiple unresolved stars, and (v) The ionizing photon flux may not be related to the star bolometric 
luminosity through standard stellar atmospheres models. For instance, it has been found that high 
accretion rates (M ~ 10 -3 — 10 -4 M y r -1 ), implies large radius a nd low effective temperatures 
associated with the accreting protostar (IHosokawa Qmukail |2009| ). High-mass protostars may 
produce little if any ionizing photons in their earliest stages. This explanation has recently been 
consider ed to account of the la ck of ionized gas towards high- mass protostellar objects in the RMS 
sample (jMottram et al.ll2011bl ). 



Table Q] presents our sample of 33 jet candidates assoc iated with high-ma ss young stars that 
fulfill the three main requisites: 23 from RMS and 10 from I Walsh et al.l (|1998l ). Columns (2) and 
(3) give the right ascension and declination of the radio source, respectively; cols. (4) and (5) the 
observed flux densities at the low and high frequencies, respectively. The low frequency corresponds 
to 4.8 or 6.7 GHz, as explained above. Column (6) give the spectral index; col. (7) the kinematic 
distance; cols. (8) and (9) the name and estimated FIR luminosity of the associated IRAS source, 
respectively; and col. (10) the reference. Interestingly, one of the selected jet candidates corresponds 



to G343. 1262— 00.0620, the most luminous sou rce known to be associated with a jet (IGaray et al 
20031 : brooks et aflbooal : lllodriguez etaDbonfJ ). The list of 33 jet candidates is an important step 



for the systematic search for jets toward HMYSOs, and the analy sis of the sub-sample presented 
in this work form the basis of the work presented in iGuzmanl ||201l( R We expect the objects listed 
in Table [1] correspond to HMYSOs in their earliest phase as radio emitting objects. 

Figure [T] presents a plot of the radio luminosity versus the bolometric luminosity for all jet 
candidates. The continuous line shows the relation between these quantities for an optically thin 
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HII region ionized by a ZAMS star (|Panagialll973l ). Sources below the dashed line fulfill the third 
selection criteria. This figure clearly illustrates that the jet candidates are underluminous in radio 
continuum emission compared to that expected for a uniform optically thin H II region. 



3. OBSERVATIONS 

The observations were made using the Australia Telescope Compact Array (atcaI) between 
2008 June and 2009 February. We combined data obtained using the 1.5B, 1.5C, and 6. OA km 
array configurations, utilizing all six antennas and covering east-west baselines from 30 m to 5.9 
km. Observations were made at four frequencies: 1.384, 2.368, 4.800, and 8.640 GHz, each with a 
bandwidth of 128 MHz, full Stokes. Throughout this work we will refer to these frequencies as 1.4, 
2.4, 4.8 and 8.6 GHz, respectively. The total integration time at each frequency was between 60 and 
240 minutes depending on each source, obtained from 10-minute scans taken over the maximum 
range of hour angles to provide good (u,v) coverage. The phase calibrators were observed for 3 
min. before and after every on-source scan in order to correct the amplitude and phase of the 
interferometer data for atmospheric and instrumental effects. The flux density was calibrated by 
observing 1934-638 (3C84) for which values of 14.95, 11.59, 5.83, and 2.84 Jy were adopted at 1.4, 
2.4, 4.8, and 8.6 GHz, respectively. 



Standard calibration and data reduction were performed using MIRIAD (jSault et al.lll995l ). 
Whenever the peak emission associated to the source was above 50 mJy beam -1 , we applied self- 
calibration on the phases. Maps w ere made by Fourier transformation of the robust-weighted 



visibilities (Robust parameter= 0, see lBriggslll995l ) . obtaining synthesized beams of typical FWHM 



of 1.6" x 1.0", 3.2" x 1.8", 6.6" x 3.5", and 10" x 5.0" at 1.4, 2.4, 4.8, and 8.6 GHz, respectively. The 
dirty maps obtained were deconvolved using the CLEAN task within MIRIAD. 

The noise levels attained in the deconvolv ed maps at 8 . 6 GH z were typically 0.1-0.2 mJy 



beam 1 , which is less than that of the RMS and lWalsh et al.l ( 19981 ) sample by factors of about 3 



and 5, respectively. Table [2] gives the phase tracking center, the on-source total integration times, 
the FWHM axes of the synthesized beam obtained at each frequency, and the deconvolved map 
noise levels for each observed field. 

4. RESULTS 

In this section we present the results of the observations towards six of the 33 jet candidates 
presented in Table [TJ The observed parameters of all the radio sources detected in 80" x 80" fields 
centered on the target source are given in Table Columns (2) and (3) give the peak position, 
cols. (4-7) the flux densities, and cols. (8-11) the deconvolved angular sizes. 



2 The Australia Telescope Compact Array is funded by the Commonwealth of Australia for operation as a National 
Facility managed by CSIRO. 
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The derived parameters, obtained from model fits to the spectra, are presented in Table [H 
Usually the spectra are well fitted by the theoretical sp ectrum from an uni form source of free- free 
emission, for which we have used the formulae given in iGarav et al.l (|1993l ). Columns (2) and (3) 
gives, respectively, the assumed distance to the source and the electron temperature used in the fit 
to the spectra. Columns (4) and (5) gives, respectively, the emission measure (EM) and angular 
size derived directly from the fit. Column (6) gives the rate of ionizing photons needed to excite 
the HII region (iVj); col. (7) the physical diameter, col. (8) the mean electron density, and col. (9) 
the type of each radio source. 

Additionally, we present in the Appendix the result of observations made toward five HMYSOs, 
selected in an early phase of this work, which are in the radio catalogs but do not fulfill the final 
jet candidate selection criteria. 

4.1. Jet candidates 



In this section we present maps of radio continuum emission maps obtained at frequencies of 
1.4, 2.4, 4.8, and 8.6 GHz, from regions of 80"x80" in size centered at the position of e ach je t 
ca ndidate observed, i ndicated by the cross, as reported in the catalogs of lUrquhart et al.l (|2007al ) 



or 



Walsh et al.l (| 19981 ). Contour levels are in a geometric progression with base between 1 and 2, 



the value chosen so each radio map displays the most important resolved and unresolved features. 
The noise and beam size for each of these maps are listed in Table [2j Squares in the 1.4 GHz maps 
indicate the peak position of the MSX sources in the field. 

In order to better determine the nature of the radio sources, we present for each field a three- 
color mid-infrared (MIR) image made using 8.0, 4.5, and 3.6 /im IRAC data from the Spitzer- 
GLIMPSE survey. The color code for these three filters is red, green, and blue, respectively. The 
brightness scales and color stretches are consistent throughout this work. We also present 24 /iin 
MIPS images obtained from the Spitzer-MIPSGAL survey for the sources G305. 7984— 00.2416 and 
G352. 5173— 00.1549. In addition, in order to support the physical interpretation developed for 
G337. 4032— 00.4037, we show a CO J = 3 — >• 2 line spectrum obtained from the APEX (Atacama 
Pathfinder Experiment) telescope. 

Figures [2] to Q3] include the radio maps, and all mid-infrared supporting figures. Radio con- 
tinuum spectra of the jet candidates are presented in Fig. [15j Figure [TBI displays spectra of other 
sources detected in the fields. We also searched the literature for the presence of water (23.2 GHz) 
and methanol (6.7 GHz) ma sers, thought to trace yo u ng stellar environment s. W e used the catalogs 



of wat er masers presented in lValdettaro et a 



(1999); and the methanol maser catalogs of 



d200lh:IUrquhart et al.l ( bopflh andlForster fc Caswell 



Pestalozzi et al 



2005) and Walsh et al 



0993 , 



199 



In the cases when a reliable position can be assigned to maser emission (via interferometer obser- 
vations) we indicate it in the radio maps. 



G305.7984 — 00.2416- This jet candidate, taken from the I Walsh et al.l survey, is associated 
with IRAS 13134-6342, although it is located ~ 24" north of the peak position of the IRAS 
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object. Figure [2] shows that there are four radio objects within the 80" x 80" region: three compact 
components (labeled A, B, and C) and an extended source, labeled D, located about 35" south of 
component A. Components A, C and D were detected at all frequencies, whereas component B was 
only detected at the two higher frequencies. Components B and C are unresolved at all frequencies. 

In order to better determine the nature and evolutionary stage of these radio sources we present 
in the top panel of Figure [3] the three color MIR image, and in the bottom panel an image of the 
Spitzer-MIPSGAL 24 //m emission. Superimposed in these two images are contours of the 4.8 
GHz emission. A comparison between the radio and MIR images shows that radio component A 
is associated with an extended MIR feature which is prominent at 8.0 /im, components B and C 
are associated with compact MIR objects GLIMPSE G305. 7996-00.2438 and G305. 7991-00.2447, 
respectively, and component D is associated with an extended envelope of diffuse emission most 
promi nent at 8.0 /im, tho ught to arise from a photon dominated (or photo-dissociation) region 
(PDR, iHeitsch et al"1l2007h . 



The jet candidate corresponds to component A, which at 8.6 GHz exhibits a cometary-like 
structure with a head toward the north and a tail trailing to the south. Its radio continuum 
spectrum (see Figure [T5|) is well fitted by that of a uniform density H II region with an emission 
measure of 2.3 x 10 6 pc cm -6 and an angular size of 1.4". The ionizing photon flux needed to excite 
this H II region is 9.6 x 10 45 photons s . We conclude that component A corresponds to a cometary 
UCH II excited by a B1-B0 .5 ZAMS star. Methanol maser activity has been detected toward this 
field (IPestalozzi et al.ll2005l ). although it was not possible to determine its precise location. 



The radio continuum spectra of components B, C and D are shown in Figure [TBI The spectra 
of source B shows that its flux density is rising with frequency, suggesting that it corresponds to 
an optically thick HII region. A fit with a uniform density HII region model gives an emission 
measure of 1.3 x 10 8 pc cm~ 6 and an angular size of 0.066". For component C we derived an 
emission measure of 1.0 x 10 6 pc cm -6 and an angular size of 0.9". The ionizing photon rate 
needed to excite this HII region is iVj ~ 2 x 10 45 s _1 , which can be provided by a Bl ZAMS star. 
For component D we derived an emission measure of 2.0 x 10 5 pc cm -6 and an angular size of 8.1". 
The ionizing photon rate needed to excite this extended and diffuse H II region is iVj ~ 3 x 10 46 s _1 , 
which can be provided by a B0. 5 ZAMS star. 

The G305. 7984— 00.2416 region may be representative of what is likely to see within the ma- 
ternities of high-mass stars: a cluster of young stars in different evolutionary stages. We find that 
the target radio source in the field is not a jet but an UC H II. The weakest radio source detected in 
the region (Component B) is particularly striking because, as shown in Fig. [3l it is associated with 
a green object i n the three color Spitz e r image, possibly indica ting shock activity. Extended green 
objects (EGOs, IChambers et al.l l2009i : ICvganowski et al.l 120081 ) are thought to probe shocked gas 
regions associated with HMYSO s through vibra tional H2 and CO bands. However, hydrogen Bra 
contribution cannot be ruled out (|Qiu et al.ll201ll ). Even though a consistent scale have been applied 
to make all three-color IRAC images presented in this work, precaution is needed in interpreting 
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the images. The community has not adopted an standardized procedure to make these images, 
casti ng doubts on the physical significance of some of the "green fuzzies" (e.g. iDe Buizer &: Vacca 
20ld ). 



The position of Component B is consistent with the peak of the emission in the MSX band E 
(21.34 /im) image (marked with a square in Fig. [2]), and it also seems to be the brightest source in 
the 24 fim MIPS image (Fig. [3]), despite the saturation of the central pixels. This evidence suggests 
that component B is actually the most embedded, and possibly the youngest, object in the region. 
Using the GLIMPSE flux densities and assuming that the MSX E ba nd flux arises ent i rely f rom this 
source, the YSO spectral energy distribution (SED) model fitter of iRobitaille et al.l (|2007l ) gives a 
total luminosity of 3.3 ±0.9 x 10 3 L & , equivalent to that of a B2 ZAMS star. The ratio between the 
ionizing photon flux derived from the radio SED fitting (JVj ~ 1.2 x 10 45 s _1 ) and that inferred from 
the total luminosity (iV* ~ 6.3 x 10 44 s _1 ) is 1.9, implying that this radio source is slightly under- 
luminous in radio. Given its characteristics, particularly the association to green fuzzy emission, 
we suggest that component B could represent the transition phase as an early B-type star grows 
onto an O-type star. The excess emission at 4.5 /xm may arise from shocks produced in the infalling 
gas. 

G311. 4298-00. 5612.- This jet candidate is associated with the RMS/MSX source G317.4298-00.5612 
(also IRAS 14477—5947). The kinematical distance of the object is 15 kpc, with no ambiguity 



( V LSR = 27.4 km s in the CS(2— >1) line, iBronfman et al. 
curve (G = 220 km s" 1 ). Fr om the bolometric flux given by | 
a total luminosity of L ~ 3.2 x 10 5 Lq, equivalent to that of a ZAMS sta r with an spectral typ e 



1996) and assuming a flat rotation 



Mottram et al.1 (|2011al ) we estimate 



between 06 and 05.5. Water maser has been detected toward this source bv lUrquhart et al.l ((2009). 



Figure H] shows that there are two radio objects within the 80" x 80" region: a compact 
component (labeled A), seen only at 4.8 and 8.6 GHz, and an extended component (labeled B) 
mainly seen at 1.2 and 2.4 GHz and completely resolved out at 8.6 GHz. 

Component A corresponds to the jet candidate. At 1.4 and 2.4 GHz this object is immersed 
within the more extended component so we could not measure its flux densities. In the frequency 
range between 4.8 and 8.6 GHz, the flux density rises with frequency. Since the source is unresolved, 
from these data it is not possible to discern whether this object is a jet or an optically thick H II 
region. If an H II region, then a fit to the radio continuum spectrum with an uniform density model 
(see Fig. [T5j) indicates an emission measure of 1.9 x 10 8 pc cm -6 and an angular size of 0.27". The 
rate of ionizing photons required to excite this region of ionized gas is 7.2 x 10 47 s -1 , which can be 
provided by a 09.5 ZAMS star. 

The radio continuum spectrum of component B (see Figure [To]) suggests it corresponds to an 
extended optically thin HII region. A fit with a theoretical spectra of an homogeneous constant 
density HII region gives an emission measure of 3.5 x 10 5 pc cm -6 and an angular size of 7.6". The 
rate of ionizing photons required to excite the H II region is 9.3 x 10 47 s" 1 , which can be provided 
by a 09.5 ZAMS star. 
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Figure [5] shows the three color mid-infrared IRAC image toward this region. Component A is 
associated with a bright MIR source with colors similar to a PDR associated to an UCH II region. 

The sum of the luminosities of the ZAMS stars exciting components A and B is 8 x 10 4 L@, 
about 4 times smaller than the bolometric luminosity of the whole region. This discrepancy may 
be solved if component A is not an HII region but an ionized jet, in which core the exciting source 
could have a larger luminosity than that implied from the observed radio flux density. Further 
observations are needed to confirm this interpretation. 

G3374032- 00.4037.- This jet candidate is associated with the RMS source G337.4032-00.4037 
(also IRAS 16351—4722). Figure [6] shows that there are two radio sources within the 80" x 80" 
region: a bright compact component (labeled A) and a weak component, labeled B, located about 
13" southeast of component A. Component A, which corresponds to the jet candidate, is unre- 
solved at all frequencies, whereas component B is resolved out at 4.8 and 8.6 GHz. Methanol maser 
activity has been associated to co mponent A bv lWalsh et al.l (|l998T ). whereas no water maser was 
detected bv lUrquhart et al.l (|2009i ) at 0.25 Jy sigma level. 



Figure [7] shows the three-color IRAC image of the MIR emission detected toward this region, 
and in red contours the 2.4 GHz data. This Figure shows that component A is associated with 
the central, bright, greenish-like o bject seen in the Spi tzer image. Assuming that this star forming 
region is at a d i stance of 3.2 kpc (IFaundez et al.ll2004l ). then from the bolometric flux reported by 
Mottram et al.l (|2011al . their component "B"), we obtain a bolometric luminosity of 2.7 x 10 4 L , 
equivalent to that of BO ZAMS star. 

The radio continuum spectrum of component A (see Figure fT5"1) shows that the flux density 
steadily increases with freque ncy. A power la w fit gives an spectral in dex of 0.9 = fc 0.15 , consistent 
with a pressure co nfined iet (lRevnoldslll986l ^. Using equation (19) of iReynoldsl (|1986l ). assuming 
typical values (e.g. iGuzman et al.ll2010l ) of 45° for the inclination, 0.2 rad for the jet aperture, a 
turnover frequency of 20 GHz and a wind velocity of 300 km s —1 , we obtain a mass loss rate of 
M « 2 x 1O~ 5 M yr- 1 . 

Our interpretation of this source as a jet is supported by the association to a very high velocity 
molecular outflow. Figure [8] shows the CO J = 3 — s> 2 emission detected toward G337. 4032— 00.4037 
using the Atacama Pathfinder Experim ent (APEX). A description of the instrument and its per- 
fo rmance is given by lGiisten et al.l (120061 ) and a detailed account of these observations can be found 
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Guzmanl (|201ll ) . The presence of high- velocity gas is evident from Fig. El where the line emission 
can be detected in a ~ 120 km s _1 range. The ambient gas velocity respect to the local standard 
of rest (LSR) is -40.7 km s _1 . 

Radio continuum spectra similar to the one exhibited by this source have us ually been inter- 
preted as arising from hyper compact HII regions (HCH IIRs, 
dependence of the density with radius (e.g.,|Panagia fc Fellilll975[ 



Franco et al 



Avalos et al 



2000) with power-law 



20061 ) . For an spher- 



ical HII region in which the density follows n e oc r 13 , the flux density should exhibit a power-law 
dependence with frequency, S v oc v a , with a = (2/3 — 3.2)/(/3 — 0.5). There is no hint of an spectrum 
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turn over at high frequencies, indicating that the emission does not reach optically thin conditions 
in the whole observed frequency range. Assuming ionization equilibrium, with JY* ~ 2.6 x 10 47 
s _1 (appropriate for a ZAMS star with 2.7 x 10 4 Lq), and that the power-law index of the density 
profile is 2.5 (determined from the above relationship), the observed spectrum and source size can 
be well reproduced if the region has an inner radius of 0.0032 pc, an outer radius of 0.044 pc and 
an electron density of 3.6 x 10 5 cnr 3 at the inner radius. We note, however, that the jet and the 
HCHIIR models have many free parameters, and that these solutions are not unique. 

Figure[7]also shows that component B is associated with the bright, central part of an extended, 
diffuse MIR object, most prominent at 8.0 /im and likely to correspon d to a PDR around an 



HII region. From the bolometric flux reported bv lMottram et al.l (|2011al ). we derive a bolometric 
luminosity for component B of 8.4xl0 4 Lq. The flux density expected at 1.4 GHz from an optically 
thin HII region excited by a star with that luminosity is ~ 0.5 Jy, considerably larger than the 
measured 5 mJy This implies, as our observations indicate, that the flux density is resolved out 
by the interferometer even at 1.4 GHz. 

G345. 0061+01. 7944.- This jet candidate is associated with the RMS/MSX source G345. 0061+01. 7944 
(also IRAS 16533-4009). Figure M shows that there are two radio objects within the 80" x 80" 
region: a compact component (labeled B) and an extended source, labeled A, located about 22" 
southwest of component B. Component B, which corresponds to the jet candidate, is unresolved at 
all frequencies whereas component A is partially resolved out at 4.8 GHz and highly resolved out 
at 8.6 GHz. 

Figure [TOl presents the three-color image of the MIR emission toward this region, together with 
red contours displaying the 2.4 GHz data. This Figure shows that component B is associated with a 
bright, compact MIR source with enhanced 4.5 (im band emission (green color) whereas component 
A is associated an extended so urce of diffus e emis sion most prominent at 8.0 /im. Methanol maser 



activity has been detected by I Walsh et al.l (|1998l ) toward two positions: one associated to radio 



component B, and another located approximately 20" to the northwest. These positions correspond 
to the posit ion of the "green" source s seen in Fig. [TUJ Water maser activity has been detected in 



the field by lForster h Caswelll (jl999l ). but not directly associated to any radio or IR source. 



The radio continuum spectrum of component B (see Figure fl~5"|) suggests it corresponds to an 
optically thick H II region. A fit with a theoretical spectra of an homogeneous constant density H II 
region, gives an emission measure of 4.6 x 10 8 pc cm -6 and an angular size of 0.73". The rate of 
ionizing photons required to excite the H II region is 1.7 x 10 47 s -1 , which can be p rovided by a BO 



ZAMS star. Component B is associated with one of the three IR sources reported bv lMottram et al 



( 2011al ) within the region (their component "C"). From its bolometric flux, of 1.19 x 10~ 7 erg s _1 
cm -2 , and assuming it is at the near distance of 1.7 kpc, we derive that this IR source has a 
luminosity of = 1.1 x 10 4 Lq , equivalent to that of a B0. 5 ZAMS star. The ionizing photon 
flux provided by such a star (JY* ~ 1.7 x 10 46 s _1 ) is about an order of magnitude smaller than 
that implied from the modeling of the radio spectra, being this source overluminous in radio rather 
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than under luminous. We note however that there is considerable uncertainty in the determination 
of the bolometr ic flux: taking a value closer to the upper limit of the uncertainty range given by 



Mottram et all (1.94 x 10" 
photoionized region. 



erg s" 



_1 cm 2 ), makes the radio and IR luminosities consistent with a 



The radio continuum spectrum of component A indicates that it corresponds to an optically 
thin region of ionized gas. A fit with a theoretical spectra of an homogeneous constant density 
HII region, gives an emission measure of 5.0 x 10 5 pc cm -6 and an angular size of 8.0". The rate 
of ionizing photons required to excite this HII region is 2.2 x 10 46 s -1 , which can be pro vided 
by a BO. 5 ZAMS star. Two of the three MIR sources identified by IMottram et al.l (|2011al . their 
components A and B) are associated with this radio component, but none is located at the peak 
of the radio emission. 



G352. 5173-00. 1549- This jet candidate, taken from the lWalsh et al.l ([ 19981 ) survey, is associ- 
ated with IRAS 17238—3516. Figure [TTl shows that there are two radio sources within the 80" x 80" 
region: a bright compact component (labeled A) detected at the four frequencies, and a weaker, 
presumably extended component (labeled B) only seen at 1.4 and 2.4 GHz and resolved out at the 
higher frequencies. 

Component A, which corresponds to the jet candidate, can be well modeled by a homogeneous 
HII region with T e ~ 15000 K, EM ~ 2.8 x 10 8 pc cuT 6 and an angular size of 0.42" (Figure 
[T5"]) . The rate of ionizing photons required to excite this HII region is 3.2 x 10 47 s _1 , which can 
be provided by a O9.5-B0 ZAMS star. The IRAC image shown in the upper panel of Figure 
PT21 shows that component A is associated with a deeply embedded object seen conspicuously in 
MIPS at 24 fj,m (lower panel), but undetected in the GLIMPSE images, indicative of very high 
absorption. Note that the object located 3" to the west of component A, seen at 4.5 and 3.6 (j,m. 
(GLIMPSE 352.5170-00.1540), and in 2MASS bands (2MASS17271105-3519310), is a foreground 
object. Assuming that the emission detected in the IRAS bands (IRAS 17238—3516) corresponds 
to reprocessed emission from the exciting source of component A and that it is located at its 
kinematical (near) distance of 5.2 kpc, we derive a total luminosity of ~ 10 5 L@, equivalent to 
that of an 07 ZAMS star. The ionizing photon rate expected from a star of that luminosity 
(N+ ~ 4.2 x 10 48 s" 1 ) is about 13 times greater than that derived from the fit to the radio spectra, 
indicating that this source is intrinsically under luminous in radio. Water maser emission has been 
detected towards the central ra dio source (IForster fc Caswelllll999i ). Methanol maser activity has 
been detected toward this field (jPestalozzi et al.ll2005l ). but its position is not well constrained by 
the observations. 



We note that the flux at 8.64 GHz reported by IWalsh et al 



is ~ 6 times smaller 

than the flux measured b y us at 8.6 GHz. Moderate flux variations in hypercompact H II regions 
have been detected (e.g. iGalvan-Madrid et al. 2008 ). but according to recent models, it is very 
unlikely to explain the larg e flux increment implied by our observations in a ~ 10 year basis 
( Galvan-Madrid et al]|2~01lh . Either the reason of the difference being observational or physical, 
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we do not discuss it further in the present work. 



Component B corresponds probably to a more developed H II region, which is consistent with 
the radio emission at the lower frequencies. Figure [121 shows that this source is associated with a 
faint 24 fim counterpart. 



G009. 9937-00. 0299.- This jet candidate, taken from lWalsh et all (|l998T ), is associated with the 
MSX source G009.9983-00.0334 and with IRAS 18048-2019. It is the only radio source within the 
80" x 80" region shown in Fig. [I3j Its radio continuum spectrum (see Fig. PT5j) can be well modeled 
with the theoretical spectra of an homogeneous constant density H II region. Assuming an electron 
temperature of 8000 K, we derive an emission measure of 5.9 x 10 5 pc cm -6 and an angular size of 
1.5". The rate of ionizing photons required to excite the H II region is 7.9 x 10 45 s _1 . Assuming that 



this object is located a t 5 kpc — the near distance derived f rom the observed vlsr 
( Bronfman et al.lll99a ) — its total luminosity is 1.7xl0 4 L & (jLackingtonl 12011 
of B0 and a B0. 5 ZAMS stars, and for which AT* = 7.9 x 10 46 s" 1 (interpolated from |Panagialll973l ) 



49 km s 1 
in between those 



The G009. 9937— 00.0299 object is then underluminous in radio by a factor of ~ 10. However, we 
do not see neither the spectral nor morphological features to characterize this radio source as a jet. 
Moreover, the three-color MIR image presented in Fig. [14] resembles that of G317. 4298— 00.5612, 
exhibiting a single central point source with color characteristic of photoionized gas. We conclude 
that this radio source corresponds to an UC HII region. 



Walsh et al.l (|1998l ) reported methanol maser emission located approximately 45" toward the 



southwest of the radio source and probably not directly r elated to it. Water maser activity has also 
been detected toward this source (Valdettaro et al. 2001). 



DISCUSSION 



5.1. Nature of the jet candidates 



In this section we summarize the nature of 7 out of the 33 jet candidates listed in Table 1. We 
find that two candidates — G337.4032-00.4037 A and G345.4938+01.4677 — can be identified, 
from their radio continuum spectral characteristics, as ionized jets. Even though the angular 
resolution of the observations was insufficient to resolve the jet morphology, there is good evidence, 
such as the presence of collimated lobes of radio emission and energetic molecular outflows to 
consider these two objects as bona-fide collimated ionized winds. 

To determine the physical p arameters of the jet f ound toward G345. 4938+01. 4677 we assumed 
a pressure confined jet model (jGuzman et al.ll2010l ). This model reproduces well the observed 
spectral indexes of G345.4938+01.4677 and G337.4032-00.4037 A. Despite of the fact that pressur e 
has been rejected as the collimation mechanism of jets associated to low-mass stars (ICabritl 120071 ). 
it is still a plausible mechanism in the high-mass case. Jets associated to HMYSOs are found 
near the center of high-mass star forming cores, which have enough density and temperature to 
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confine a jet with the characteristics described above. Furthermore, ram pressure from infalling 
gas is also relevant in these cores. Infalling motions have been detected towa rd all cores associated 
with luminous HMY SOs with jets: G343.1262 -00.0620. iGarav et al.l l20ol: G345.4938+01.4677, 



Guzman et all 1201 ll : and G337.4032-00.4037 A ( Guzman 



20111 ) 



Three of the jet candidates, G317.4298-00.5612 A, G345. 0061+01. 7944 B, and G352. 5173-00.1549 
A, turned out to correspond to small {D < 0.03 pc) regions of ionized gas with large emission mea- 
sures (EM > 10 8 pc cm~ 6 ), and thu s can be classified as hypercompact HII regions, according to 
the definition of iSewilo et al.l (|201ll ). HCHII regions are thought to trace the earliest stages of 
photoionized nebulae associated to s ingle or binary hig h-mass star system system. To date, about 
a dozen of these regions are known (jHoare et al.l 120071 ) . Finally, the remaining two jet candidates 
observed, G009.9937-00.0299 and G305. 7984-00.2416 A, correspond to UCHII regions (EM> 10 6 
pc cm" 6 , size < 0.1 pc). An important issue concerning the jet list candidates is whether the 
radio and FIR emission arise from the same object. At least in one case, G305. 7984— 00.2416 A, we 
found that the cataloged RMS radio source is not associated with the luminous object seen in IRAS 
bands, and that the much fainter component B is likely to be the youngest and more embedded 
source of this high-mass star maternity. 

C onsidering the 7 je t candidates observed as part of this search, plus G343. 1262—00.0620, 
which Gar ay et al. (|2003l ) showed that corresponds to a collimated ionized jet, then there are 8 of 
the 33 jet candidates for which their nature have been determined. Of these, 37.5% corresponds to 
collimated ionized winds, 37.5% corresponds to HC HII regions and 25% corresponds to UCHII 
regions. Due to the homogeneity of the selection criteria used in building the jet candidate list 
(Table [1]) we expect the same proportion of objects in the whole list. Thus, we expect to find about 
9 new ionized jets among the 25 candidates not yet observed. 



5.2. Lifetime of jets in HMYSOs 



Despite of the limited size of the sample of objects observed in this work, we can analyze the 
statistical incidence of jets from our observations and draw some conclusions about the lifetime of 
the jet phenomena in the formation process of high-mass stars. For the statistical analysis we will 
only consider the jet candidates chosen from the RMS survey. This selection allows us to make 
comparison with the results of a recently reported study of HM YSOs and compact H II regions 
(CHllRs) identified in the RMS survey (jMottram et al.l l2011bl ). In particular, they estimated 
that the lifetime of the CHIIR phase, for all range of luminosities, is ~ 3 x 10 5 yr. Of the 239 



radio s ources we initially considered from the entire RMS radio sample reported bv lUrquhart et al 
(|2007al ). we find that 92 are located at angular distances smaller than 25" from an IRAS point 



source with a luminosity L > 2 x 10 4 L@. We assume that the s e 92 so urces form an unbiased 
sample of the CHII region population analyzed in lMottram et al.l (|2011bl ). 



Of these 92 sources we find that 23 fulfill the two additional criteria required to be considered 
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jet candidates: positive spectral index and radio underluminosity by a factor of at least 10. These 
are the objects with entry (1) in col. (10) of Table [TJ Since we can not disentangle the FIR 
luminosity from the jet candidates G301. 1364-00.2249 A and B, and G317.8908-00.0578 A and B, 
we count them as if they were single objects. At the high luminosity of the sources we are sampling, 
completeness corrections of the RMS survey are small: about 0.8 and 0.97 for so urces of 20, 000 



and 3 0, 000 Lq, respectively, and essentially complete for more luminous HMYSOs (jMottram et al 



201 lbh . 



We observed four jet candidates drawn from the RMS subsample and found that two of them 
are bona-fide jets (G337.4032-00.4037 A and G345.4938+01.4677). Using this detection rate, we 
then expect to find about 11 jets within the 23 candidates. The jet incidence in the population of 
CH II regions is therefore ^11/92, w hich in combination with the CH II region average lifetime of 



3 x 10 yr (jMottram et al.ll2011bl ) implies that the lifetime of the jet phase is roughly 4 x 10 yr. 



This lifetime is comparable to the Kelvin- Helmholtz timescale of a HMYSO of 10 4 Lq. Afterwards, 
the central protostar rapidly contracts onto the main sequence producing an HCH II region. This 
could explain why few jets are observed associated with HMYSOs. 

Similarly, for the HCH II regions we derive lifetimes of the order of 4 x 10 4 yr. Multiple radio 
sources in the field is an issue to consider within the RMS sample, but we found that all companion 
sources correspond to more diffuse and evolved H II regions. 



5.3. The role of jets in high-mass star formation 



Our systematic search have shown that jets occur in the formation process of high-mass stars, 
at least for HMYSOs with luminosities up to 1 x 10 5 Lq. This result strongly suggests a disk- 
mediated accretion scenario for the formation of high-mass stars. The duration of the jet phase is 
short, ~ 4 x 10 yr, explaining why they are not commonly observed. 

This age estimate allows us to assess the role of the collimated ionized jets in driving the 
larger scale molecular outflows o bserved toward HMYSOs. The mean mass loss rate estimated 

is ~ 5 x 10~ 6 Mr. 



from the values of known jets (see iGaray et al 



2003 



2006 



Rodriguez et al. 



2008 



Guzman et al. 



Assuming that the jet veloc i ties are ~ 500 km s" 1 (|Angladalll99i 



2011 



Marti et al. 



1998; 



Curiel et al 



Guzman et al.ll2010l ). then the momentum delivered by the jets is 



~ 100 M & km s . On the other hand, th e mean masses and gas velocities of the massive molecular 
outflows reported bv lBeuther et al.l ( 2002 ) are 25 Mq and 28 km s , respectively, implying a mean 
momentum of ~ 700 Mq k m s -1 . Thus, even with the assumption that all of the jet momentum is 
transferred to the outflow ([Richer et al.ll2000l ). the jet momentum as estimated above falls short, 
by a factor of ~ 7, to drive the molecular o utflow. This factor is similar to the one reported for 
G345. 4938+01. 4677 by iGuzman et al.1 (1201 ll ). There are at least two possible answers to explain 
this discrepancy: i) the jet is not entirely ionized, but their ionization fraction is about a 10% or 
below; ii) as appears in the radio images, the ejection of mass from the jet is episodic and highly 
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discontinuous, most of the momentum being delivered during short bursts, and what we observe 
now correspond to a more "quiescent" phase of the jet. 



6. SUMMARY 

We are currently undertaking a systematic search of collimated jets towards high- mass young 
stellar objects with the goal of determining the pertinence of this phenomenon in the formation 
process of massive stars. In this paper we present the first results of multi-frequency radio contin- 
uum observations, made using ATCA, towards 6 HMYSOs candidates to harbor jets. Our main 
results and conclusions are summarized as follows. 

1. We compiled a list of 33 HMYSOs likely to contain jets using a selection criteria based on their 
radio and infrared properties. The objects have large bolometric luminosities (L > 2 x 10 Lq) 
and are underluminous at radio wavelengths compared to what is expected from the total 
luminosity. 

2. We obtained radio continuum observations at 1.4, 2.4, 4.8, and 8.6 GHz using ATCA, towards 
six HMYSOs in the list. We find that one object corresponds to an ionized jet, three to 
hypercompact H II regions and two to ultracompact H II regions. 

3. The jets discovered as part of this search, G337.4032-00 . 4037 a nd G345.4938+01.4677 - the 



latter reported and analyzed previously by I Guzman et al.l (|2010l ) - are associated with the two 



most luminous (7x 10 4 and 1.0 x 10 5 Lq) HMYSOs known to harbor this type of objects. This 
indicates that the phenomena of collimated ionized winds appears in the formation process 
of stars at least up to masses of ~ 20 M and provides strong evidence for a disk-mediated 
accretion scenario for the formation of high-mass stars. Prom the rate of occurrence of jets 
in our sample, we estimate that the jet phase in high-mass protostars lasts for ~4x 10 yr. 

4. The estimated momentum that the observed ionized jets deliver during its main lifetime is 
about one tenth of the average momentum of the molecular outflows associated to HMYSOs. 
Whether this is a problem of the models or an actual inadequacy of the jets to explain the 
dynamics of molecular outflows remains to be resolved by future observations. 

5. An important part of the observed candidates were found to correspond to hypercompact H II 
regions, thought to be an early phase in the development of the UCH IIR. These hypercompact 
regions probably appears right after the jet phase. We estimated their lifetime in a similar 
manner as done with the jets, obtaining also ~ 40, 000 yr. 

6. As part of this search, we performed observations toward 5 other HMYSOs that did not 
satisfied the criteria to be jet candidates. In the Appendix we present these observations and 
physical interpretation of the data. 
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Other sources of interest 



In this sect i on we rep ort the resul t s of t he observations on five radio sources, chosen from the 
Urquhart et al.l (|2007al ) or I Walsh et al.l (|1998l ) catalogs, but that do not fulfilled all the requirements 
to be considered jet candidates. These sources were selected in an early phase of our study, and 
they are most likely HMYSOs. These sources were observed together with the jet candidates, and 
the same reduction and calibration steps described in § [3] were applied. 

Table Ogives the phase tracking center, the on-source total integration times, the FWHM axes 
of the synthesized beam obtained at each frequency, and the deconvolved map noise levels for each 
observed field. 

Emission was detected toward four of the five objects: three of them correspond to ionized 
regions associated to a HMYSO, and one to an extragalactic source. Maps of the radio continuum 
emission at the frequencies of 1.4, 2.4, 4.8, and 8.6 GHz, in regions of 80"x80" in size, are shown 
in Figures [17] through [20] The position of the target radio object, as reported in either catalog 
of is indicated in the maps with a cross. We also show in Figure [21] a mid-infrared 3 color image 
taken toward G333. 1306—00.4275, that will be used to support its physical interpretation. Spectra 
of the four detected sources are shown in Fig. [22] The derived parameters, obtained from model 
fits to the spectra, are presented in Table [7] The columns display the same parameters as Tabled] 
described in § 0] 

G263. 7759- 0042813.- This radio source is associated with the RMS /MS X object G263.7759- 00.4281 
(also IRAS 08448—4343 ), which is known to b e associated with an H2O maser ( Urquhart et all2009 ) 
and a H2 bipolar flow (jGiannini et al.ll2005l ). The lumi nosity derived from th e IRAS bands and 
assuming this source is located at a distance of 1.6 kpc (jUrquhart et al.ll2007bl ). is ~ 3 x 10 3 L , 
therefore failing our luminosity selection criterion. It is the only radio source within the 80" x 80" 
region shown in Fig. [TTJ Its radio continuum spectrum (see Figure [22]) can be well modeled with 
the theoretical spectra of an homogeneous constant density H II region. Assuming an electron tem- 
perature of 8000 K, we derive an emission measure of 4.1 x 10 6 pc cm -6 and an angular size of 
0.53". The rate of ionizing photons required to excite the HII region is 7.0 x 10 44 s _1 . 



Using MSX, IRAS and TIMMI2 data, iMottram et alJ (j2011al ) derived a bolometric flux of 
Fho\ ~ 5.77 x 10 -8 erg cm -2 s _1 . At a distance of 1.6 kpc, we derive a bolometric luminosity of 



L ~ 4.6 x 10 3 L Q , equivalent to that of B1-B2 zero-age mam seq 



flux from such a star is 1.4 x 10 45 s 1 (interpolated from iPanagial 1 19731 ) . similar to that required to 



uence (ZAMS) star. The ionizing 
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excite the region of ionized gas. We conclude that this radio object corresponds to a small (0.004 
pc in diameter), optically thin HII region excited by a single B1-B2 ZAMS star. 

G268. 6162-00. 7389.- This radio source is associated with the RMS/MSX source G268.6162-00.7389 
(also IRAS 09014-4736). It is the only radio source within the 80" x 80" re gion shown in Fig. [TSJ 
No water masers were detected toward this source by lUrquhart et al.l (J2009j) at the 0.25 Jy sigma 
level. The l uminosity derived from the IRAS bands and assuming this source is located at a distance 



of 1.7 kpc ([Beck et al.lll99ll ). is ~ 5 x 1O 3 L0, therefore failing our luminosity selection criterion. 



The radio source is unresolved and its radio continuum spectrum (see Figure [22]) is well fitted by 
that of a uniform density HII region with an emission measure of 7.0 x 10 6 pc cm -6 , temperature 
of 15000 K and an angular size of 0.84". The ionizing photon flux needed to excite the HII region 
is 3.4 x 10 45 photons s _1 . 



Mottram et al.l (|2011al ) reported for this object a bolometric flux of 7.1 x 10 W m . At a 



distance of 1.7 kpc, the implied bolometric luminosity is 6.4 x 10 3 L & , equivalent to that of a ZAMS 
star with a spectra l type between B0. 5 and Bl. The ionizing flux from such a star is 3.4 x 10 45 s _1 
(interpolated from iPanagial Il973l ). equal to that required to excite the region of ionized gas. We 
conclude that this radio object corresponds to a small (0.007 pc in diameter), optically thin HII 
region excited by a single B0.5-B1 ZAMS star. 

G327. 1307+00. 5259.- This radio source is in the field of the RMS source MSX-G327.1192+00.5103, 
although it is located at an angular distance of ~ V to the north of the peak MSX emission and 
the closest IRAS point source (IRAS 15437—5343). No infrared counterpart has been detected at 
MIR wavelengths, neither in Spitzer (GLIMPSE or MIPSGAL) nor MSX, which makes it unlikely 
that this object corresponds to an embedded HMYSO. 

It is the only radio source detected within the 80" x 80" region shown in Fig. [191 Its radio 
continuum spectrum (see Figure l22j) shows a decrease in the flux density at the higher frequencies. 
Being unresolved at all wavebands, the drop is not produced by the resolving power of the inter- 
ferometer. We suggest that the radio emission from this source has a non-thermal origin. From 
the available evidence, we further suggest that this is an extragalactic source, possible a radio-loud 
galaxy. Figure [22] sh ows that the s pectrum can be well fitted with a simple model of a gigahertz 
peaked galaxy ([Snellen et al.lll998l ). the spectral indeces in the optically thick and thin regimes 
being 0.91 and —0.69, respectively. The lack of infrared emission could be explained by the source 
being intrinsically faint at these wavelengths compared to the radio emission, and/or due to ab- 
sorption in the Galactic plane. Similar objects without MIR co unterpart had already been detected 
in extragalactic surveys (jNorris et al.ll2006l : iHuynh et al.ll2010l ). 



G333. 1306— 00.4275 - This radio object was selected from the lUrquhart et al.l catalog due to 
its steep spectral index (> 1) between 4.8 and 8.6 GHz. It is located ~ 25" east of the closest MSX 
source G333. 1306—00.4275 (also IRAS 16172—5028), therefore failing our selection criteria. There 
are two radio sources within the 80" x 80" region shown in Fig. [20j a compact central component 
(the target object; labeled B) and a complex, extended component, labeled A, associated with the 
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RMS source. 



The radio continuum spectrum of component B (see Figure [22]) is we ll fitted with a pow er-law 
spectrum with a spectral index of ~ 1, simil ar to those of HC H II region ([Franco et al.ll2000l ). The 
bolometric flux is 15 .4xl0 -8 erg s -1 cm -2 ( Mottram et al.ll2011al ). which at a distance of 3.5 kpc 
( Faundez et al.l 120041 ) implies that its total bolometric luminosity is 5.9 x 10 4 L & . 



Figure [211 shows a three-color image of the MIR emission towards G333. 1306—00.4275 made 
using the 3.6, 4.5, and 8.0 /im data obtained from the Spitzer/GLIMPSE survey, and superimposed 
in red contours the radio emission detected in 4.8 GHz. This figure clearly shows that component 
B is associated with a green extended object, thought to be related to regions of shocked gas. Com- 
ponent A is associated with a bright, extended MIR source, specially bright at 8.0 //m, associated 
with PAH emission in a PDR. The radio spectrum is flat (Fig. [22]) and probably correspond to a 
more evolved HII region. 



G345. 3768+01. 3926.- This source was chosen from the catalog of IWalsh et all (11998T). It 



is located more than 1' away from the nearest IRAS source (IRAS 16561—4006). IWalsh et al 



reported a peak flux of 5 mJy beam -1 at 8.64 GHz and a flux densitjjfl of 12.3 mJy No sources 
were detected towards the radio position in GLIMPSE or MSX, which makes this object an unlikely 
YSO candidate. We do not detect emission towards this radio object in any of the four frequencies 
at the ~ 0.1 mJy beam -1 level, indicating that this is a highly variable radio source. We propose an 
extra galactic origin. Such extreme variability in radio-loud galaxies is rare, but has been observed 



(e.g., iBarvainis et al.ll2005l ). 



The integrated fluxes were published online in VizieR jOchsenbein et aljbooch 



Table 1. Jet candidates 



Source Name 


a 


5 


F„_low a F 8 . 6G Hz S.I. D 


IRAS 


Liras Ref. 




(J2000) 


(J2000) 


(mJy) (mJy) (kpc) 




(10 4 L ) 



G240 


3160- 


fOO 


0714 


07 h 


44 m 52?04 


-24°07 


42'.'4 


<18 


G274 


0649 


-01 


1460 


09 


24 


42.13 


-52 


02 


00 


8 


28.9 


G289 


9446 


-00 


8909 


11 


01 


09.00 


-60 


56 


56 


3 


< 0.7 


G293 


9633 


-00 


9776 


11 


32 


36.14 


-62 


28 


08 


3 


24.8 


G298 


2234 


-00 


3393 


12 


10 


01.16 


-62 


49 


53 


9 


1350. 


G300 


9674+01 


1499 


12 


34 


53.23 


-61 


39 


40 


1 


102. 


G301 


1364 


-00 


2249A 


12 


35 


35.13 


-63 


02 


31 


7 


83.8 


G301 


1364 


-00 


2249B 


12 


35 


35.19 


-63 


02 


24 





126. 


G305 


7984 


-00 


2416 


13 


16 


42.62 


-62 


58 


21 


2 


<10 


G308 


9176+00 


1231 


13 


43 


01.72 


-62 


08 


56 


1 


247. 


G309 


9196+00 


4791 


13 


50 


41.89 


-61 


35 


11 


5 


377. 


G311 


1359 


-00 


2372 


14 


02 


09.93 


-61 


58 


37 


9 


3.0 


G317 


4298 


-00 


5612 


14 


51 


37.60 


-60 


00 


19 


4 


6.6 


G317 


8908 


-00 


0578A 


14 


53 


06.19 


-59 


20 


56 


7 


< 0.7 


G326 


4477 


-00 


7485 


15 


49 


18.67 


-55 


16 


52 


5 


4.7 


G328 


5759 


-00 


5285B 


15 


59 


38.15 


-53 


45 


27 


9 


326 


G331 


4181 


-00 


3546 


16 


12 


50.24 


-51 


43 


28 


6 


80.7 


G333 


0162+00 


7615 


16 


15 


18.70 


-49 


48 


52 


8 


46.9 


G336 


9842 


-00 


1835 


16 


36 


12.42 


-47 


37 


58 





18. 


G337 


4032 


-00 


4037 


16 


38 


50.45 


-47 


28 


02 


7 


90.3 


G337 


7051 


-00 


0575 


16 


38 


29.63 


-47 


00 


35 


3 


76.3 


G337 


8442 


-00 


3748 


16 


40 


26.67 


-47 


07 


13 


1 


11.1 


G340 


0708- 


fOO 


9267 


16 


43 


15.69 


-44 


35 


16 





48.8 


G340 


2768 


-00 


2104 


16 


48 


53.30 


-45 


10 


22 


3 


< 1 


G343 


1262 


-00 


0620 


16 


58 


17.21 


-42 


52 


07 


1 


<7.3 


G345 


0061- 


fOl 


7944 


16 


56 


47.59 


-40 


14 


25 


8 


127. 


G345 


4938- 


fOl 


4677 


16 


59 


41.61 


-40 


03 


43 


4 


4.8 


G352 


5173 


-00 


1549 


17 


27 


11.32 


-35 


19 


32 


8 


<14 



12.4 


u 


7.3 


07427 


-2400 


7.7 


(2) 


36.1 


0.44 


6.9 n 


09230 


-5148 


80. 


(1) 


12.4 


> 2 


10.3 


10591 


-6040 


21.5 


(1) 


28.9 


0.3 


11.1 


11303 


-6211 


10.1 


(1) 


2200. 


0.96 


11.3 


12073 


-6233 


463. 


(1) 


138. 


0.6 


4.4 


12320 


-6122 


22.6 


(1) 


228. 


2. 


4.4 


12326 


-6245 


32. 


(1) 


179. 


0.68 


4.4 


12326 


-6245 


32. 


(1) 


5. 


u 


3. 


13134 


-6242 


3.4 


(2) 


374. 


0.81 


5.3 


13395 


-6153 


25.6 


(1) 


384. 


0.072 


5.4 


13471 


-6120 


27. 


(2) 


3.5 


0.3 


14.3 


13585 


-6144 


11.7 


(1) 


8.2 


0.42 


15.0 


14477 


-5947 


63.4 


(1) 


3.6 


> 2 


15.1 


14492 


-5908 


19.7 


(1) 


5.9 


0.39 


4.3 


15454 


-5507 


2.1 


(1) 


381 


0.27 


11.4 


15557 


-5337 


368. 


(1) 


83.9 


0.07 


4.1 


16090 


-5135 


5.7 


(1) 


88.8 


1.2 


3.3 


16115 


-4941 


6.1 


(1) 


34.3 


1.3 


10.8 


16325 


-4731 


38. 


(1) 


130. 


0.71 


3.2 


16351 


-4722 


10.6 


(1) 


171. 


1.6 


12.2 


16348 


-4654 


53.1 


(1) 


24.4 


1.5 


3.1 


16367 


-4701 


3.9 


(1) 


93.8 


1.1 


4.9 


16396 


-4429 


7.2 


(1) 


4.7 


> 2 


3.6 


16452 


-4504 


8.0 


(1) 


4.62 


u 


2.9 


16547 


-4247 


6.3 


(2) 


209. 


0.98 


1.7 n 


16533 


-4009 


5.4 


(1) 


12.5 


1.9 


1.7 


16562 


-3959 


7.0 


(1) 


7.5 


u 


5.6 n 


17238 


-3516 


9.9 


(2) 



Table 1 — Continued 



Source Name 


a 


5 


Fu-low a Fg.6GHz S.I. D 


IRAS 


L/_ras R-ef. 




(J2000) 


(J2000) 


(mJy) (mJy) (kpc) 




(10 4 L Q ) 



G000. 3138-00. 2000 
G009.9937-00.0299 
G010.8403-02.5913 
G024.4673+00.4910 
G025.6469+01.0534 



17 47 09.66 -28 46 

18 07 52.84 -20 18 
18 19 12.10 -20 47 
18 34 08.12 -07 18 
18 34 20.91 -05 59 



27.7 


8.3 


14.5 


29.3 


<10 


2.32 


30.7 


<21 


4.8 


18.2 


<28 


13 


39.3 


<7 


1.4 



2.2 


8.0 n 


17439 


u 


5.0 n 


18048 


u 


1.9 n 


18162 


u 


5.9 


18314 


11 


3.1 


18316 



2845 


48. 


(2) 


2019 


2.1 


(2) 


2048 


2.5 


(2) 


0720 


44. 


(2) 


0602 


2.9 


(2) 



Here i^-low represents 4.8 or 6.7 GHz depending whether the candidates comes from the RMS or from lWalsh et al.l 



(|1998l ) catalog, respectively. 

"Distance ambiguity not resolved. Near distance adopted. 



References. — (1 nUrquhart et al.l |2007al) . (2): IWalsh et all |l998l ) 



Table 2. Observational parameters of the observed jet candidates 













Synthesized Beam 






Noise (mjy beam x ) 




Source 


Phase Tracking Center 


Time 


1.4GHz 


2.4GHz 


4.8GHz 


8.6GHz 


1.4GHz 


2.4GHz 


4.8GHz 


8.6GHz 




a(J2000) 


<5(J2000) 


(mins.) 


(") 


(") 


(") 


(") 


(mjy) 


(mjy) 


(mjy) 


(mJy) 


G305. 7984-00. 2416 


13 16 42.62 


-62 58 21.2 


90 


6.6 x 6.3 


4.4 x 3.9 


2.2 x 1.7 


1.1 x 0.9 


0.2 


0.2 


0.07 


0.08 


G317.4298-00.5612 


14 51 37.60 


-60 00 19.8 


110 


9.9 x 5.4 


6.8 x 3.8 


2.8 x 1.7 


2.1 x 1.1 


0.2 


0.2 


0.1 


0.08 


G337.4032-00.4037 


16 38 50.45 


-47 28 02.7 


240 


8.1 x 6.0 


5.4 x 3.9 


2.2 x 1.7 


1.1 x 0.9 


0.4 


0.2 


0.08 


0.09 


G345. 0061+01. 7944 


16 56 47.59 


-40 14 25.8 


180 


8.7 x 5.9 


5.6 x 4.0 


2.6 x 2.0 


1.4 x 1.0 


0.3 


0.3 


0.1 


0.2 


G352. 5173-00.1549 


17 27 11.32 


-35 19 32.8 


120 


10.1 x 5.6 


6.4 x 3.7 


3.1 x 2.0 


1.7 x 1.0 


0.2 


0.2 


0.2 


0.1 


G009. 9937-00.0299 


18 07 52.84 


-20 18 29.3 


110 


20.1 x 5.1 


12.5 x 3.0 


5.9 x 1.7 


3.2 x 0.9 


0.2 


0.1 


0.06 


0.05 



Table 3. Observed parameters of radio sources 





Coordinates 




Flux density 






Deconvolved Sizes 






R.A. 


Dec. 


1.4 GHz 


2.4 GHz 4.8 


GHz 


8.6 GHz 


1.4 GHz 


2.4 GHz 


4.8 GHz 


8.6 GHz 


Source 


(J2000) 


(J2000) 


(mjy) 


(mjy) (mjy) 


(mjy) 


(") 


(") 


(") 


(") 




Jet candidates 


G305.7984-00.2416 A 


13 h 16 m 42?63 


-62°58'20'/9 


10.4 


11.2 


11.2 


9.3 


U 


1.4 


1.4 


1.4 


G317.4298-00.5612 A 


14 51 37.63 


-60 00 20.2 






11.0 


19.7 






U 


0.36 


G337.4032-00.4037 A 


16 38 50.47 


-47 28 03.1 


27.8 


56.7 


110.3 


139.5 


1.6 


1.2 


0.5 


0.8 


G345.0061+01.7944 B 


16 59 37.76 


-40 12 03.8 


12.2 


44.4 


152.5 


271.8 


1.3 


U 


0.7 


0.7 


G352.5173-00.1549 A 


17 27 11.32 


-35 19 32.2 


4.4 


14.1 


40.5 


66.0 


U 


u 


0.9 


0.4 


G009.9937-00.0299 


18 07 52.82 


-20 18 28.9 


3.4 


3.5 


3.2 


2.3 


U 


u 


U 


U 






Additional sources detected 


in the fields 










G305.7984-00.2416 B 


13 16 43.23 


-62 58 32.9 






0.7 


1.0 






u 


U 


G305.7984-00.2416 C 


13 16 43.41 


-62 58 28.7 


1.2 


2.2 


1.8 


2.0 


U 


u 


u 


0.9 


G305.7984-00.2416 D a 


13 16 44.50 


-62 58 54.7 


32.6 


27.6 


24.2 


35.0 


8.0 


10 


10 


10 


G317.4298-00.5612 B a 


14 51 38.46 


-60 00 24.3 


41.6 


48.9 


5.7 




7.9 


7.1 


7 




G337.4032-00.4037 B a 


16 38 51.13 


-47 28 14.8 


5.0 


5.1 


6.5 


4.7 


U 


U 


2.6 


~ 2 


G345. 0061+01. 7944 A a 


16 56 45.97 


-40 14 38.1 


91.4 


99.1 


83.0 


72 


7.1 


8.1 


~ 8 


~ 8 


G352. 5173-00.1549 B a 


17 27 11.91 


-35 19 35.5 


2.8 


1.6 






U 


U 







a The peak positions for these sources was obtained from the 1.4 GHz images. 



-24- 



Table 4. Derived parameters of the radio sources 





Dist. 


T e 


EM 


e 


N» 


D 


(«e> 


Type 


Source 


(kpc) 


(10 4 K) 


(pc cm~ 6 ) 


(") 


(lO^s- 1 ) 


(PC) 


(cm -3 ) 






Jet candidates 


G305. 7984-00. 2416 A 


3.0 


0.8 


2.30 x 10 6 


1.40 


9.56 


0.020 


1.1 x 10 4 


UC HII 


G317.4298-00.5612 A 


15. 


0.8 


1.90 x 10 8 


0.27 


718 


0.020 


9.8 x 10 4 


HC HII 


G337. 4032-00. 4037 A 


3.2 


0.8 


> 2 x 10 s 


< 0.64 




< 0.01 


> 1.5 x 10 5 


jet 


G345.0061+01.7944 B 


1.7 


1.5 


4.60 x 10 8 


0.73 


170 


0.006 


2.8 x 10 5 


HC HII 


G352. 5173-00. 1549 A 


5.2 


1.5 


2.80 x 10 8 


0.42 


318 


0.011 


1.0 x 10 5 


HC HII 


G009.9937-00.0299 


5.0 


0.8 


5.94 x 10 5 


1.49 


7.89 


0.036 


4.1 x 10 3 


UC HII 






Additional sources 


detected 


in the fields 






G305. 7984-00. 2416 B 


3.0 


1.0 


1.3 x 10 8 


0.066 


1.23 


0.001 


3.67 x 10 5 


HC HII 


G305. 7984-00. 2416 C 


3.0 


0.8 


1.0 x 10 6 


0.901 


1.74 


0.013 


8.74 x 10 3 


UC HII 


G305. 7984-00. 2416 D 


3.0 


0.8 


2.0 x 10 5 


8.027 


27.6 


0.117 


1.31 x 10 3 


C HII 


G317.4298-00.5612 B 


15.0 


0.8 


3.5 x 10 5 


7.030 


926. 


0.511 


8.30 x 10 2 


C HII 


G337.4032-00.4037B 


3.2 


0.8 


3.0 x 10 5 


2.657 


5.16 


0.041 


2.70 x 10 3 


UC HII 


G345.0061+01.7944 A 


1.7 


0.8 


5.0 x 10 5 


8.027 


22.2 


0.066 


2.75 x 10 3 


UC HII 


G352. 5173-00. 1549 B 


5.2 


0.8 


5.0 x 10 5 


1.329 


5.68 


0.033 


3.86 x 10 3 


UC HII 



Table 5. Observational parameters of other sources of interest 













Synthesized Beam 






Noise (mjy beam x ) 




Source 


Phase Tracking Center 


Time 


1.4GHz 


2.4GHz 


4.8GHz 


8.6GHz 


1.4GHz 


2.4GHz 


4.8GHz 


8.6GHz 




a(J2000) 


<5(J2000) 


(mins.) 


(") 


(") 


(") 


(") 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


G263.7759-00.4281 


08 h 46 m 34!85 


-43°54'29?8 


60 


14.0 x 5.0 


9.0 x 3.3 


4.6 x 1.9 


2.5 x 1.0 


0.1 


0.1 


0.07 


0.08 


G268. 6162-00.7389 


09 03 09.51 


-47 48 27.3 


60 


13.0 x 4.4 


9.2 x 3.2 


4.7 x 1.9 


2.6 x 1.0 


0.3 


0.1 


0.06 


0.07 


G327. 1307+00.5259 


15 47 32.47 


-53 51 30.9 


180 


7.2 x 6.3 


4.9 x 4.2 


2.3 x 1.9 


1.2 x 1.0 


0.2 


0.1 


0.09 


0.07 


G333. 1306-00.4275 


16 21 02.95 


-50 35 12.3 


60 


12.1 x 4.4 


7.0 x 2.7 


4.9 x 1.5 


2.8 x 0.8 


3.1 


1.8 


2.1 


3.4 


G345. 3768+01.3926 


16 59 37.75 


-40 12 03.5 


170 


9.0 x 5.6 


5.7 x 3.7 


2.6 x 2.0 


1.4 x 1.0 


0.1 


0.1 


0.08 


0.08 



Table 6. Observed parameters of other radio sources of interest 





Coordinates 




Flux density 






Deconvolved Sizes 






R.A. 


Dec. 


1.4 GHz 


2.4 GHz 4.8 GHz 


8.6 GHz 


1.4 GHz 


2.4 GHz 4.1 


3 GHz 


8.6 GHz 


Source 


(J2000) 


(J2000) 


(mjy) 


(mjy) (mjy) 


(mjy) 


(") 


(") 


(") 


(") 


G263.7759-00.4281 


08 46 34.848 


-43 54 30.27 


2.0 


3.2 2.7 


2.4 


U 


U 


U 


U 


G268. 6162-00.7389 


09 03 09.499 


-47 48 27.66 


8.6 


9.1 10.0 


9.2 


U 


1.6 


1.1 


0.7 


G327.1307+00.5259 


15 47 32.472 


-53 51 31.58 


20.0 


22.5 19.5 


13.3 


u 


0.9 


u 


0.36 


G333. 1306-00.4275 B 


16 21 02.921 


-50 35 12.92 


186 


604 1100 


2000 


u 


U 


5.9 


3.3 


G345.3768+01.3926* 




















Additional sources detected in the fields 


G333. 1306-00.4275 A 


16 21 00.30 


-50 35 08.30 


1100 


1589 1348 


1664 




~ 25 







"This source was not detected in our observations. 
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Table 7. Derived parameters of the radio sources 



Source 


Dist. 

(kpc) 


(10 4 K) 


EM 

(pc cm~ 6 ) 


6 Nj 
(") (lO^s- 1 ) 


D 

(PC) 


(n e ) 
(cm- 3 ) 


Type 










Other sources 


of interest 








G263. 7759-00.4281 
G268. 6162-00.7389 
G333. 1306-00.4275 B 


1.6 
1.7 
3.5 


0.8 
1.5 
1.5 


4.05 x 10 6 
7.00 x 10 6 

> 2 x 10 8 


0.53 0.697 
0.84 3.40 
<2.5 


0.004 
0.007 
< 0.042 


3.2 x 10 4 
1.9 x 10 4 
> 7.1 x 10 4 


UC HII 
UC HII 
HC HII 



Additional sources detected in the fields 



G333. 1306-00.4275 A 3.5 0.8 1.0 x 10 6 24.998 1820 0.424 1.54 x 10 3 CHII 
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Fig. 1. — Radio luminosity at 8.6 GHz versus Liras luminosity (Eq. [T]) for all jet candidates in 
Table [TJ The continuous line indicates the expected radio luminosity of an optically thin H II 
region ionized by a ZAMS star fPanagia 1973 ). The dashed line illustrates graphically the third 
selection criterion. Sources observed as part of the jet search described in this work are encircled. 
Squares mark objects G343. 1262-00.0620 (also IRAS 16547-4247), G345.4938+01.4677 (IRAS 
16562-3959), and G337.4032-00.4037. 
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Fig. 2. — ATCA maps of the radio continuum emission toward G305. 7984— 00.2416. The cross 
marks the position of the radio source reported by Walsh et al. (1998). Beams are shown in the 
lower left corner of each panel. Top left: 1.4 GHz map. The squares mark the peak position of two 
MSX 21 fim sources in the field. Top right: 2.4 GHz map. Bottom left: 4.8 GHz map. Labeled are 
the four radio sources detected in this field. Bottom right: 8.6 GHz map. Contour levels are —5, 
5, 8, 13, 20, 30, 46, and 70 times 0.16 mJy beam -1 for the 1.4 and 2.4 GHz maps, and 0.07 mJy 
beam -1 for the 4.8 and 8.6 GHz maps. 
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Fig. 3. — Comparison between the radio and IR images for G305. 7984— 00.2416. The white and 
black contours in each panel show 4.8 GHz data. Top panel: 3 color IRAC image using 8.0, 4.5, 
and 3.6 fim data from Spiteer-GLIMPSE for red, green, and blue. Bottom panel: 24 fim MIPS 
image from Spiteer-MIPSGAL. Saturated pixels are shown in magenta. 
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Fig. 4. — ATCA maps of the radio continuum emission toward G317. 4298— 00.5612. The cross 
marks the position of the radio source reported by Urquhart et al. (2007a). Beams are shown in 
the lower left corner of each panel. Top left: 1.4 GHz map. The square marks the peak position of 
the MSX source detected at 21 ^m. Top right: 2.4 GHz map. Bottom left: 4.8 GHz map. Labeled 
are the two radio sources detected in this field. Bottom right: 8.6 GHz map. Contour levels are 
-5, 5, 7, 9, 13, 19, 27, 40 times 0.5, 0.5, 0.2, and 0.4 mJy beam" 1 for the 1.4, 2.4, 4.8 and 8.6 GHz 
images, respectively. 
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Fig. 5. — Comparison between mid-infrared and radio emission data detected toward 
G317. 4298— 00.5612. The color image correspond to a three color IRAC image using 8.0, 4.5, 
and 3.6 /im data from 5piteer-GLIMPSE for red, green, and blue. Blue contours show 4.8 GHz 
data. 
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Fig. 6. — ATCA maps of the radio continuum emission toward G337. 4032— 00.4037. The cross 
marks the position of the radio source reported bv Urauhar t et al. (2007a). Beams are shown in 
the lower left corner of each panel. Top left: 1.4 GHz map. The squares show the peak position 
of the three MSX 21 (iva. sources in the field. The star marks the position of the methanol masers 
detected by Walsh et al. (1998). Top right: 2.4 GHz map. Labeled are the two radio sources 
detected in the field. Bottom left: 4.8 GHz map. Bottom right: 8.6 GHz map. Contour levels are 
-5, 5, 14, 28, 50, 81, 129, and 200 times 0.25 mJy beam" 1 for the 1.4 and 2.4 GHz images, and 
-5, 5, 29, 65, 118, 199, 319, and 500 times 0.15 mJy beam" 1 for the 4.8 and 8.6 GHz images. 
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Fig. 7.— Three color MIR image toward G337.4032-00.4037 made using Spitzer-GLIMPSE images 
at 8.0, 4.5, and 3.6 /j,m for red, green, and blue, respectively. Red contours correspond to the 2.4 
GHz radio data. 
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Fig. 8.— CO(3->2) spectrum observed toward G337.4032-00.4037 using APEX. The v LSR of the 
ambient gas associated to the source is —40.7 km s _1 . High- velocity molecular gas emission is 
detected between —85 and 41 km s _1 , spanning 126 km s _1 . 
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Fig. 9. — ATCA maps of the radio continuum emission towards G345. 0061+01. 7944. The cross 
marks the position of the radio source reported by Urauhart et al. (2007a) . Beams are shown in 
the lower left corner of each panel. Top left: 1.4 GHz map. The square marks the peak position 
of the MSX source detected at 21 //m, the triangle the position of the H2O maser reported by 
Forster Caswell (1999), and the stars the position of methanol masers detected by Walsh et al. 
(1998). Top right: 2.4 GHz map. Bottom left: 4.8 GHz map. Labeled are the two radio components 
detected in the field. Bottom right: 8.6 GHz map. Contour levels are —5, 5, 11, 19, 31, 50, 78, 120 
times 0.4 mJy beam -1 for the 1.4 and 2.4 GHz images and -5, 5, 14, 28, 50, 81, 129, 200 times 0.4 
mJy beam -1 for the 2.4 and 8.6 GHz images. 
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Fig. 10.— Three color Spitzer MIR image of G345.0061+01.7944, made using 8.0, 4.5, and 3.6 
//m data from 5piteer-GLIMPSE for red, green, and blue, respectively. Superimposed red contours 
show the 2.4 GHz radio emission. 
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Fig. 11. — ATCA maps of the radio continuum emission towards G352. 5173— 00.1549. Beams are 
shown in the lower left corner of each panel. The cross marks the position of the radio source 
reported by Walsh et al. (1998 ). Top left: 1.4 GHz map. Contour levels are —5, 5, 8, 13, and 
18 x 0.20 mJy beam -1 . Labeled are the two radio components detected in the field. Squares 
mark the position of MSX sources. The triangle marks the position of the H2O maser detected by 
Forster fc Caswell (1999). Top right: 2.4 GHz map. Contour levels are -5, 5, 8, 13, 18, 26, 36, 
and 50 x 0.20 mJy beam -1 . Bottom left: 4.8 GHz map. Contour levels are —5, 5, 8, 13, 18, 26, 
36, and 50 x 0.60 mJy beam -1 . Bottom right: 8.6 GHz map. Contour levels are —3, 3, 9, 18, 29, 
44, 64, and 90 x 0.60 mJy beam -1 . 
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Fig. 12. — Mid-infrared images toward G352. 5173— 00.1549, with black contours corresponding to 
1.4 GHz emission as seen in Fig. II 1L Top Panel: 3 color IRAC image using 8.0, 4.5, and 3.6 fim data 
from 5piteer-GLIMPSE for red, green, and blue, respectively. The peak of the 1.4 GHz emission 
is consistent with the peak detected at higher frequencies, and to radio component A. The source 
seen ~ 3" to the west is not embedded, and is unrelated with the radio emission. Radio component 
B appears to be related with the 8.0 ^m diffuse emission seen at south-east from the peak. Bottom 
Panel: MIPS 24 /im emission. Component A is consistent with the peak at MIR wavelengths, and 
component B also shows MIR counterpart. 
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Fig. 13. — ATCA maps of the radio continuum emission toward G009. 9937— 00.0299. The cross 
marks the position of the radio source reported bv Walsh et al. (1998). Beams are shown in the 
lower left corner of each panel. Top left: 1.4 GHz map. The square marks the peak position of 
the associated MSX source G009. 9983— 00.0334 and the star the position of the methanol masers 
detected bv Walsh et al. (1998 ). Top right: 2.4 GHz map. Bottom left: 4.8 GHz map. Bottom 
right: 8.6 GHz map. Contour levels are —3, 3, 4, 5, 7, and 10 times 0.30 mJy beam -1 for the 1.4 
and 2.4 GHz maps, and 0.23 mJy beam -1 for the 4.8 and 8.6 GHz maps. 
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Fig. 14. — Comparison between the mid-infrared and radio emission detected toward 
G009. 9937— 00.0299. The color image correspond to a three color IRAC image using 8.0, 4.5, 
and 3.6 /xm data from iSpiteer-GLIMPSE for red, green, and blue. Blue contours show 4.8 GHz 
data. 
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Fig. 15. — Radio continuum spectra of the jet candidates observed in this work. The source name 
is given in the upper left of each panel. The dashed line corresponds to the best-fit obtained with 
an homogeneous HII region model. The derived emission measure (in pc cm" 6 ) and solid angle (in 
arcsec 2 ) are indicated in the lower right of each panel. Dotted lines in the G337. 4032— 00.4037 A 
panel show the best fit obtained with ionized jet models (Reynolds 1986). 
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Fig. 16. — Left column: Spectral energy distributions of G305. 7984— 00.2416 B, C and D. 
Right column: Spectral energy distributions of G317.4298-00.5612 A, G337.4032-00.4037 B, and 
G345. 0061+01. 7944 A. The dashed line shows the best-fit obtained with an homogeneous HII re- 
gion, with the emission measure and solid angle indicated in each panel in units of pc cm -6 and 
arcsec 2 , respectively. 
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Fig. 17. — ATCA maps of the radio continuum emission towards G263. 7759— 00.4281. The cross 
marks the position of the radio source reported bv Urauhar t et al. (2007a). Beams are shown in 
the lower left corner of each panel. Top left: 1.4 GHz map. The squares marks the peak position 
of the MSX sources. Top right: 2.4 GHz map. Bottom left: 4.8 GHz map. Bottom right: 8.6 GHz 
map. Contour levels are —3, 3, 5, 7, and 10 times 0.15 mJy beam -1 for the 1.4 GHz image and 
—3, 3, 5, 7, 10, and 14 times 0.15 mJy beam -1 for the 2.4, 4.8, and 8.6 GHz images. 
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Fig. 18. — ATCA maps of the radio continuum emission towards G268. 6162— 00.7389. The cross 
marks the position of the radio source reported bv Urauhar t et al. (2007a). Beams are shown in 
the lower left corner of each panel. Top left: 1.4 GHz map. A square marks the position of the 
MSX source. Top right: 2.4 GHz map. Bottom left: 4.8 GHz map. Bottom right: 8.6 GHz map. 
Contour levels are —5, 5, 7, 9, 11, and 15 times 0.40 mJy beam -1 for the 1.4, 4.8, and 8.6 GHz 
images, and —5, 5, 7, 9, 11, 15 and 20 times 0.40 mJy beam -1 for the 2.4 GHz image. 
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Fig. 19. — ATCA maps of the radio continuum emission toward G327. 1192+00. 5103. The cross 
marks the position of the radio source reported by Urquhart et al. (2007a ) . Beams are shown in 
the lower left corner of each panel. Top left: 1.4 GHz map. Top right: 2.4 GHz map. Bottom left: 
4.8 GHz map. Bottom right: 8.6 GHz map. Contour levels are —5, 5, 8, 13, 18, 26, 36, 50 times 
0.35 mJy beam -1 for the 1.4 and 2.4 GHz images, and 0.2 mJy beam -1 for the 4.8 and 8.6 GHz 
images. 
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Fig. 20. — ATCA maps of the radio continuum emission toward G333. 1306—00.4275. The cross 
marks the position of the radio source reported by Urquhar t et al. (2007a). Beams are shown in 
the lower left corner of each panel. Top left: 1.4 GHz map. Squares mark the positions of MSX 
sources. Top right: 2.4 GHz map. Bottom left: 4.8 GHz map. The two main radio components 
are marked with letters A and B. Bottom right: 8.6 GHz map. Contour levels are —7, 7, 17, 31, 
49, 74, 107, and 150 times 3 mJy beam . 
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Fig. 21.— Three color IRAC image of G333. 1306-00.4275 made using 8.0, 4.5, and 3.6 [an data 
from Spitzer-GhlMPSK for red, green, and blue, respectively. Intercalated contours of the 4.8 GHz 
data (Fig. [20]) are shown in red. 
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Fig. 22.— Le/f column: Spectral energy distributions of G263.7759-00.4281, G268.6162-00.7389, 
and G327. 1192+00.5103. In the first and second panel the dashed line shows the best-fit obtained 
with an homogeneous HII region, with the emission measure and solid angle indicated in each 
panel, in units of pc cm~ 6 and arcsec 2 , respectively. In the third panel the dot-dashed line is the 
best fit obtained with a gigahertz peaked source model f Snel len et al. 1998 s ). Right column: SEDs 
of G333. 1306—00.4275 A and B. In the second panel (component B), the dashed line shows the 
best fit obtained with a rising power-law. 



